Gaped graphene bilayer: disorder and magnetic field effects 
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Double layer graphene is a gapless semiconductor which develops a finite gap when the layers are placed 
at different electrostatic potentials. We study, within the tight-biding approximation, the electronic properties 
of the gaped graphene bilayer in the presence of disorder, perpendicular magnetic field, and transverse electric 
field. We show that the gap is rather stable in the presence of diagonal disorder. We compute the cyclotron 
effective mass in the semi-classical approximation, valid at low magnetic fields. Landau level formation is 
clearly seen in zigzag and armchair ribbons of the gaped bilayer at intermediate magnetic fields. 



Graphene, the two-dimensional allotropic form of carbon, has recently deserved considerable attention. The realization that 
graphene could be obtained and studied experimentally revealed a plethora of unusual properties that may be useful in the design 
of new electronic devices [1-3]. The peculiar band structure of graphene, where at the corners of the Brillouin zone (BZ) the 
excitations are massless Dirac fermions, is responsible for many of the unconventional properties ]4|]. 

The possibility of creating stacks of graphene layers with the accuracy of a single atomic layer, providing an extra dimension 
to be explored, is another advantage of graphene for electronic applications. In particular, bilayer graphene has shown to have 
unusual electronic properties, though unexpectedly dissimilar to those exhibited by its single layer parent. The new type of 
integer quantum Hall effect observed in bilayer graphene fUEI, which is induced by chiral parabolic bands, is an example of its 
uniqueness. Recent advances on the experimental side made it possible to produce bilayers where the two layers are effectively 
at different electrostatic potentials |7l l8[] - introducing the concept of biased bilayer. This asymmetry between layers opens a 
gap between valence and conduction bands, which, in the unbiased case (both layers at the same potential), touch in a parabolic 
way at zero energy. The gap is a function of the potential between layers, and may be controlled by changing the electric field 
across the bilayer. This switching functionality opens the door for potential applications of bilayer graphene on atomic-scale 
electronic devices |9|]. 

In this paper we study the biased bilayer system starting from the simplest tight-binding model. By computing the density of 
states (DOS) we analyze how the gap structure evolves when diagonal disorder is present. In the presence of a perpendicular 
magnetic field we obtain the cyclotron effective mass (semi-classical approximation), and study the Landau level formation in 
zigzag and armchair ribbons of bilayer graphene. The paper is organized as follows: in Sec. [TT] we present the tight-binding 
Hamiltonian we use and discuss basic aspects of the electronic structure; the effect of disorder in the gap structure is discussed 
in Sec. [TTTJ in Sec. [IV] we analyze the effects of a perpendicular magnetic field; and Sec. [Vl contains our conclusions. 



In single layer graphene the carbon atoms form a honeycomb lattice which has two atoms per unit cell. These two atoms 
belong to different triangular lattices labeled A and B. The bilayer is made of two layers, which we label 1 and 2, arranged in 
the Bernal stacking (A1-B2). The unit cell has then four atoms, one per layer per sublattice. Each carbon atom contributes with 
a 7T electron, and thus the system is naturally at half-filling. 

In the tight-binding approximation the relevant energy scales are the in-plane hopping energy, t « 2.7 eV, and the interlayer 
hopping energy, t± w 0.35 eV. When a perpendicular magnetic field B = B e z is applied to the system t± is unaffected but t 
acquires a phase IT^l such that, t — ► te 1 ^ Ir A dr , where e is the electron charge, 8 is a vector connecting nearest-neighbor 
sites, and A is the vector potential. The tight-binging Hamiltonian describing non-interacting it electrons in bilayer graphene 
then reads: 



I. INTRODUCTION 



II. MODEL AND BASIC ELECTRONIC STRUCTURE 
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FIG. 1: Band structure of single layer graphene (a), unbiased bilayer graphene (b), and biased bilayer (c) with V = 2t±, along two directions 
in the BZ near the Dirac (K) points. The inter-layer hopping was set to t±/t = 0.2. 

where a\ . , [oi,(m,n)] an ^ i m n \ \pi,{m,n)\ creates (annihilates) an electron on atom Ai and Bi of layer i at cell (m, n), 
respectively. Spin indices have been neglected for simplicity. The Peierls phase in Eq. (fl]i is written assuming a Landau gauge, 
A = {—y, 0, 0)B, where <fi/(f)o — BAq /<fio is the magnetic flux through a plaquette in units of the flux quantum <fia = h/e (Aq 
being the area of the graphene unit cell). The last term in Eq. (HJ stands for the electrostatic potential difference, V, between the 
two layers, with nAi,( m ,n) ar, d n Bi,(m,n) number operators. 

The first (or second) line of Eq. (Q]l describes tt electrons on single layer graphene. The resultant dispersion relation is given 
by, 

4A 2 = 3 + 2cos(ak x ) + 4cos(afc a /2)cos (ak v V3/2\ , (2) 

where a is the in-plane C-C distance for sites of the same sublattice. The vertex of the valence and conduction bands touch at 
the corners of the BZ (Dirac points), as shown in Fig.[T](a). A Dirac linear dispersion follows at low energies, e(q) — ±vpHq, 
where = taK~ 1 \/3/2 is the Fermi velocity which substitutes the speed of light, and q is the wave vector relatively to the 
Dirac points. Taking into account the full Hamiltonian given in Eq. (Q]i we obtain four bands with dispersion energies given by, 

E^(V)=±^ei+tl/2+Vyi± y£i/4 + (tl+V^ei , (3) 

with ek as in Eq. @. For the unbiased case Eq. ([3]) reduces to E± — ±(e k + ± t±/2, from which it becomes 

clear that the spectrum is as in Fig.[T](b): two low energy gapless bands touching in a parabolic way at the Dirac points, and 
two high energy bands with a gap of t±. As the undoped bilayer is at half-filling, the Fermi energy (Ep) occurs at exactly at 
the points the bands Ey~ touch. Thus, low energy quasi-particles are massive Dirac fermions with effective mass given by 
to* = 2t±h 2 /3a 2 t 2 ll^ TTlll . The biased bilayer, however, shows a true gap between valence and conduction bands, as shown in 
Fig- [11(c), and the low energy bands have a "Mexican hat" like dispersion. For V <C t, which is the case for the biased bilayer, 
the gap behaves as A g = [t\V 2 /{t\ + V 2 )] 1 ^ 2 , and is therefore fully controlled by the bias V. 



III. EFFECT OF DIAGONAL DISORDER IN THE GAPED STRUCTURE 



Recent experiments using angle-resolved photoemission spectroscopy (ARPES) have clearly shown that the spectrum of 
the biased bilayer is well described by Eq. ©18], with the measured energy-momentum dispersion in good agreement with 
Fig.Q](c). However, the experimental data also shows the presence of a finite spectral weight inside the gap, which may be 
attributed to disorder. The sources of disorder in the biased bilayer are not yet understood. Inhomogeneities in the substrate, 
below the bottom layer, as well as in the dopant coverages, above the top layer, are plausible sources. Here we consider the biased 
bilayer in the presence of diagonal disorder, adding the term Hdi s — J2 m n »=i 2 £ i,i m - n ) \ n M,{m,n) + n m,{m,n)] to Eq. ©. We 
use the box distribution function to generate the random on-site energies, e G [-W/2, W/2], and vary its width W in order to 
induce spectral weight inside the unperturbed gap. 

To study how diagonal disorder affects the gap, and in particular how W compares with the other energy scales of the system, 
we have computed the DOS with the recursion method 012I1 . We define zero temperature retarded Green's functions in the 
standard way, 

G a r y(t) = -i (o| K r (t), 4 r ,(o)} |o> e(t), G b r b r ;(t) = -i (o| {b itT (t), &l )r ,(o)} |o> e(t), (4) 

where r = (m, n) specifies the cell. The recursion method gives an approximate value for the retarded Green's function in the 
thermodynamic limit, and therefore for the disordered DOS, by simulating large lattices. The disorder averaged DOS is then 
easily obtained by averaging over disorder realizations. 
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FIG. 2: (Color online) (a) - Non-disordered DOS: exact (thick lines) and recursion method result (normal line), (b) - DOS for several values 
of disorder W. (c) - The same as in (b) closer to the gap region. 



Figure |2] resumes the results we have obtained for the effect of diagonal disorder in the DOS, with emphasis in the gap 
behavior. The simulated lattices had 2d 2 sites, with d = 1280, and we averaged over 100 disorder realizations. The inter-layer 
hopping was set to t±/t = 0.2. The bias was fixed to V/t = 0.02, as experimentally it is the smallest energy scale, V <C tj_, 
and thus the gap is essentially the bias, A g w V. In order to check the performance of the recursion method in describing such a 
tiny region of the spectrum (A g <C t± <C t), we show in Fig. 12(a) the exact non-disordered DOS [13] along with the recursion 
method result. The non-disordered DOS of the unbiased bilayer (dashed line) falls linearly when \E\ — > 0, being finite at zero 
energy II 1 lL 1 1 3H . When the layers are made inequivalent by the bias a gap opens and ID-like divergences associated with the 
"Mexican hat" like dispersion region show up 11411 . The recursion method reproduces this behavior fairly well, in spite of a 
somewhat underestimated gap (doted vertical lines signal the exact the gap limits). In Fig. [2(b) we show the DOS for several 
disorder values. The disorder parameter W equals the relevant energy scales of the problem: bias, V; inter-plane hopping, t±; 
and in-plane hopping t. It is readily seen that the disordered DOS starts to depart from the non-disordered one only for W ~ t. 
Figure|2(c) is our main result of this section, and shows the disordered DOS in the gap region. Diagonal disorder closes the gap 
when its width is of the order of the largest energy scale in the system, W ~ t. 



IV. MAGNETIC FIELD EFFECTS 



We have computed the cyclotron effective mass m*(n), whose semi-classical expression is given by m*(n) = 
{fx 2 /2ir)dA(E) /dE\E=E F (n), where A{E) is the fc-space area enclosed by the orbit of energy E and n is the carrier den- 
sity at the Fermi energy Ep. Its dependence on the carrier density n may be found through measurements of Shubnikov de 
Haas oscillations. Here we consider only the experimentally relevant case V <C t± <C t, and we further assume that the in- 
equality Ep <C t holds. In that case the Dirac linear dispersion can be used for et in Eq. ©, and analytical expressions for 
the cyclotron effective mass may be derived. When E F is varied we have to distinguish three different cases (see Fig. |2: for 
|£f| in the "Mexican hat" region, |£f| < V/2, there are two types of quasi-particles with Fermi wave vectors q£ and (rel- 
atively to Dirac points) and opposite cyclotron effective masses; when |.Ef| is between the top of the "Mexican hat" and the 
bottom of the high energy band, V/2 < E F < {t\ + V 2 /4)^ 2 , only quasi-particles exist; q F quasi-particles show up again 
for Ep > (t 2 ^ + V 2 /4) 1 / 2 with positive cyclotron effective mass. The Fermi wave vectors are functions of the Fermi energy, 
V 2 /4 ± [Ep(V 2 + t\) — t\V 2 /4] 1 / 2 } 1 / 2 , and the respective cyclotron effective masses are given by, 
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with h = a 2 t 2 n. Figure [3] shows the result for m+ [Eq. ©J as a function of the carrier density; n < for holes and n > 
for electrons. As can be seen for V/t — 0.01 ~ 27 (meV) - full fine - there is a l/n behavior for small densities, which is 
associated with the presence of the "Mexican hat" dispersion lfl3ll . and for higher densities the cyclotron mass increases as the 
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FIG. 3: (Color online) Cyclotron effective mass mj [given 
in Eq. lO] in units of the bare electron mass, m e , as a func- 
tion of the carrier density for three different bias values. We 
set t±/t = 0.1 



FIG. 4: (Color online) Energy spectrum for a ribbon of bilayer 
graphene with zigzag (a-b) and armchair (c-d) edges and width 
N = 400 unit cells: (a-c) - V = 0; (b-d) - V = t±/W. We 
set B = 30 T and t± /t = 0.2. 



carrier density increases. It is worth mentioning that a more realistic calculation would account for the variations of the bias with 
the density, V(n), as varying one implies the variation of the other. We postpone to future work the problem of the determination 
of V(n) 111311 . Also, the bias may not be the same for equal concentration of holes and electrons |7l|g]. In that case the cyclotron 
effective mass is asymmetric with respect to doping with holes or electrons. This is clearly seen in Fig.[3j where we show m+ 
for V/t = 0.05 « 135 (meV) - holes - and V/t = 0.005 « 13.5 (meV) - electrons. 

Landau level formation in bilayer graphene was already studied in the continuum limit in Ref. @], Here we study the problem 
of a bilayer subjected to a perpendicular magnetic field using the tight-binding Hamiltonian given by Eq. (Q~|i, both at zero and 
finite bias. Again we assume that inequalities V <C t± <C t hold. Bilayer nano-ribbons with zigzag and armchair edges were 
diagonalized, the obtained spectrum is shown in Fig. [4] For each k, the momentum parallel to the edge, there are AN bands. In 
Fig- 111 however, we focus on the low energy behavior near the Dirac points. Panels (a) and (c) show the result for the unbiased 
bilayer with zigzag and armchair edges, respectively. Fourfold degenerate zero energy Landau levels are clearly seen, along with 
twofold degenerate non-zero Landau levels, in agreement with the continuum result [6]. The major difference between zigzag 
and armchair is that the former have zero energy surface states along with zero energy bulk Landau levels [13]. The result for 
the biased bilayer, V = t±/10 = 0.02, is shown in panels (b) and (d) for zigzag and armchair edges, respectively. The fourfold 
degeneracy of zero energy Landau levels is lifted. A gap V opens and twofold degenerate Landau levels at V/2 and — V/2 show 
up. In fact it can be shown that their wave functions are either localized in layer 1 or 2. For zigzag edges it becomes clear that 
surface states and bulk Landau levels of the same band live in different layers [13]. As a consequence, dispersive surface states 
appear inside the gap, as shown in Fig.|4](d). 



V. CONCLUSIONS 



Using a tight-binding approximation we have studied the electronic behavior of a graphene bilayer with layers at different 
electrostatic potential - biased bilayer. The applied bias opens a gap in the spectrum which is completely controlled by the 
applied voltage. We have shown that diagonal disorder reduces the size of the gap, which finally closes when the width of the 
disorder distribution equals the larger energy scale in the system: the in-plane hopping t. We have also studied the biased bilayer 
in the presence of a perpendicular magnetic field. We have calculated the semi-classical cyclotron effective mass as a function 
of the carrier density and bias, which is valid for low magnetic fields. When the field is sufficiently high Landau level formation 
in zigzag and armchair nano-ribbons is perfectly seen, where a gap between the lowest electron-like and the highest hole-like 
bulk Landau levels opens in the presence of a finite bias. 
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